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Abstract 

The fermions of the Standard Model are integrated out to obtain the effective Lagrangian in 
the sector violating P and CP at zero temperature. We confirm that no contributions arise for 
operators of dimension six or less and show that the leading operators are of dimension eight. To 
assert this we explicitly compute one such non- vanishing contribution, namely, that with three Z°, 
two W + and two W~ . Terms involving just gluons and W's are also considered, however, they turn 
out to vanish in the P-odd sector to eighth order. The analogous gluonic term in the CP-odd and 
P-even (C-odd) sector is non-vanishing and it is also computed. The expressions derived apply 
directly to Dirac massive neutrinos. All CP-violating results display the infrared enhancement 
already found at dimension six. 
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I. INTRODUCTION 



A full understanding of CP-violation remains a challenge and for this reason it is a 
fruitful field of research both in relation with the Standard Model of particle physics and in 
extensions thereof [l^. CP- violation enters in very different phenomena, like non vanishing 
of the electric dipole momentum of elementary particles, baryogenesis [10|, or, assuming 
CPT invariance, the puzzling T-violation. Yet, in the Standard Model, CP-violation is 
rather elusive. There is no trace of it in the QCD sector, while in the electroweak sector 
it enters through a small parameter in the CKM matrix for quarks and possibly also 



for leptons, for massive neutrinos 12]. Even in the electroweak sector manifestation of CP 
breaking requires a subtle combination, the Jarlskog determinant A, which requires order 
twelve in the quark (or leptons) masses and would vanish if two up-like or two down-like 



quarks were degenerated in mass (13J. In any case only through fermions CP can be broken 
in the Standard Model. The structure of the Standard Model action implies that integration 
of the fermions results in an effective Lagrangian of the form (we assume the unitary gauge 
throughout) 



d„-4 



^HE^f^yJ o a (x), (i.i) 

where O a (x) represents any possible operator, of mass dimension d a , constructed as a 
Lorentz and gauge invariant product of the gauge fields, their derivatives and derivatives 
of the Higgs field. g a is the operator coupling constant, with mass dimension 4 — d a . (f){x) 
denotes the Higgs field and v its vacuum expectation value. The coupling constant (which 
may vanish for some operators) has two additive contributions, one from the quark loop and 
another from the lepton loop. In the CP-odd sector, g a must contain the Jarlskog determi- 
nant. In terms of the Yukawa coupling this yields a tiny dimensionless number, A/v 12 , of the 
order of 10 -24 . This fact has occasionally been presented as an indication of an intrinsic lim- 
itation of the Standard Model to produce enough CP-breaking to account for observations, 
including the baryon asymmetry. While this might be true, qualitative arguments should 



eventually be supported by a detailed computation. Smit argued in 14] that the coupling g a 
is just a homogeneous function of the quarks (or leptons) masses of the appropriate degree. 
This implies that g a ~ A x I a , where I a has a large negative degree to compensate that of 
A. Both A and I a depend only on the fermion masses and do not involve v. On the other 
hand, the various quark masses are very different and widely different result can be obtained 



by combining them at random. Actual calculations have been carried out in [15|, [16( for op- 
erators of dimension six, which is the first possible CP-violating contribution at one-loop. 
They show that q n ~ JK/m 2 c where is m c the charm quark mass, J = 2.9(2) x 10~ 5 is the 



Jarlskog invariant 



17] and k is a dimensionless coefficient of the order of unity. Implications 



for cold electroweak baryogenesis have been considered in 18l. Il9j]. Unfortunately, these two 



references differ in that 
whereas 



15| finds such a dimension six contribution in the P-odd sector 
16] finds a contribution in the C-odd sector but none in the P-odd one. 
The purpose of this note is manyfold. First, to reduce to the simplest and more trans- 
parent terms the calculation of these couplings constants. Second, to confirm that, although 
dimension six CP-odd and P-odd operators do exist, their coupling vanish in the Standard 
Model. Third, to verify that the order six cancellation is accidental, and non vanishing con- 
tributions in the CP-odd and P-odd sector appear for the first time at dimension eight. The 
purely gluonic leading (eighth) order term is also computed since it is particularly simple. 
As it turns out, this term breaks C but not P. Lastly, to verify that the enhancement (as 
compared to the naive estimate) found at order six is displayed also at higher orders. 

II. THE METHOD 

We will integrate out the fermions in the Standard Model to extract the CP violating 
contribution of the resulting effective action. This is the one-loop approximation to the 
effective action with full one-particle irreducible bosonic lines and vertices. We work at zero 
temperature. Quarks will be explicitly considered. Leptons would not contribute to the 
CP-odd sector if neutrinos are assumed to be exactly massless. For massive Dirac neutrinos 
the contribution of the leptons will be completely analogous to the one obtained for quarks. 

The quark-sector Lagrangian of the Standard Model, in its Euclidean version and in the 

n 

unitary gauge, can be written as |20J: 



C{x) = q(x)Bq(x) = (q L ,q R ) • (2.1) 




Here qL,n carry Dirac, generation (family), ud and color indices [ud space distinguishes the 
up-like from down- like quarks in each generation). Expanding further the matrices in ud 



3 




space: 

Here m u ^ are the diagonal matrices (in generation space) with the up-like and down-like 
quarks masses, respectively. the gluon field, the Z° field, the W boson fields, C 
is the CKM matrix, finally (-D M ) u ,d = <9 M + q u ,dB IJ , where q u = 2/3, = —1/3, and is the 
weak hypercharge gauge connection. For convenience, in all cases the coupling constant has 



been included in the corresponding gauge connection. Further details can be found in 



16] . 



After integration of the quark loop, the corresponding Euclidean effective action is just 

r = -TrlogD. (2.3) 

r is the sum of the all Feynman graphs with one quark loop and any number of bosonic legs, 
gauge fields and Higgs. This sum is written as a functional which will be expressed within 
a covariant derivative expansion of these bosonic fields. 

Certainly, the effective action can be computed following the efficient method outlined in 
16] and based on 2l|, applied there to sixth order in the derivative expansion. However, 



one of our goals here is to present a derivation as transparent as possible, and closer to the 
method introduced in 22] on which the calculation of [l5| is based. To this end, we will use 
the relation 

5T = -Tr (5D D" 1 ) = - J d 4 x tr [SD^D -1 \x)] . (2.4) 

In the second equality it has been used that the variation SD, induced by the variation in the 
gauge and Higgs fields in D, contains no derivatives. The method is then to choose a suitable 
variation of these fields, compute ST in the desired sector, and subsequently seek a functional 
fulfilling such a variation. The virtues of this approach are i) the current (a;|D _1 |x) is easier 
to obtain than the TrlogD itself, ii) the condition on 5T of being a consistent variation 
provides a nontrivial check of the calculation, and iii) even if one where to compute T 
directly, the simplest way to avoid integration by parts identities (i.e., redundant operators 
in the final expression) is to obtain its functional derivative, 5T/5D = — (x|D _1 |x). This 
quantity is local and so free from ^-integration by parts identities. This issue becomes 
increasingly important as the number of derivatives increases. 
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A convenient field to use as variation in eq. ( 12. 4 j) is Z which appears just in Dr}. We 



adopt such a choice, namely, 




5p R = 6m = 0. 



(2.5) 



The Dirac operator can be written as D = PlP r Pr + P^p L P^ + P^itlPr + P^mP^ where 
Pr,l — |(1 i 75) project on the R or L spaces, respectively. D can be explicitly inverted in 
chiral blocks. In particular, for the RL block, PrD -1 Px, = Pr(P l — mp R ~ l m)~ 1 PL. Hence, 
from eq. fl23]> and fl23|) . 



5r = -Tr(P i j^( J p L -m^ i? - 1 



m 



(2.6) 



Therefore, if only the P-odd sector (i.e., that with a 75) is retained, we will have 



5r 



2 
2 



m) 



75 ^(^ + p^mp-'mp- 1 + ■ 



(2.7) 



It can be noted that each term of the expansion between parenthesis starts and ends with 
a label L. Moreover, as the chiral label propagates through the term it flips (from L to R 
and vice versa) at m but not at Dr or D^. Keeping these rules in mind we can simply write 

1. 



5T' 



-Tr 



-Tr 



75 5f{p^ + p~ x mp- x mp- x 

75 W+m)- 1 ; 



(2.* 



The second equality follows from the fact that terms with an odd number of m's are auto- 
matically discarded since they cannot start and end with a label L. 
A simp 



16 



23 



e and convenient technique to compute (x\(p+ m)~ l \x) is the method of symbols 
24} . This method is suitable for computing one- loop Feynman graphs with external 
legs at zero momentum or more generally, for expansions around zero momentum. In the 
present case it takes the form 

1 f d 4 xd 4 p 



5T- 



tr 



l5 8p(p+ii>+my l 



(2.9) 



Two remarks apply here: i) the momentum variable like Dp , does not introduce a 
flip in the chiral label, and ii) after momentum integration all P M appear only in the form 



[Dp, ] and so there are no longer differential (or pseudo differential) operators acting; only 
an ordinary function of x survives. At the same time gauge invariance is ensured. 

The co variant derivative expansion is just an expansion in powers of Dp. Only even 
orders contribute (the space-time dimension being even). Besides, in the P-odd sector the 
effective action starts at fourth order in four dimensions, since a Levi-Civita pseudo tensor 
must be present. Hence: 



5T- = i /"^tr 



2 J (2ttY 



00 / -J. \ 2n +3 • / 

^— ' \ rr + m z J rr + m z 

n=0 ' 



(2.10) 



III. SIXTH ORDER P-ODD TERMS 



Taking the appropriate values of n in eq. (I2.10p and the appropriate contributions in 
£)R,l an( j m ^ Qne can se i ec t the desired terms of the effective action. At least two W + W~ 
pairs must be present in the CP-odd sector terms, since the quark loop must visit the three 
generations [111]. 1 A term with just (W^ + W^~) 2 would count as fourth order, however no CP- 
odd term can be constructed without introducing further fields or derivatives. To sixth order 
such a P-odd and CP-odd term can be written using (W + W~) 2 ZD. (D here refers to either 
D u or Dj.) The question is whether this operator appears with a non vanishing coefficient 



in the Standard Model or not. Ref. 15] claims that it does whereas the calculation in 
concludes that it does not. Therefore we will start by reconsidering such a term within our 
present approach. 

Under a variation of Z, the contribution of the candidate term to be found in eq. f)2.10p is 
of the form SZ(W + W~) 2 D. We can set <fi = v, since we are not interested in contributions 
from Higgs, and likewise we can set Z^ and G /t to zero in D^ ,L . Moreover, we can even 
set D™ = D^ = dp in D^ ,L . The ordinary derivatives can be unambiguously replaced by 
covariant ones at the end without loss of information since no Fp V tensor can be present in 
the term considered. 

The computation is tedious but straightforward. Let us spell out the main steps in the 
calculation. We select terms with n — 1 in eq. ( 12 . 1 j) and restore the L, R labels. We keep 
only terms starting and ending with the label L and only m introduces a chiral label flip 



1 Of course, terms with a single W + W pair are allowed beyond one-loop. 
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L -H- R. No flip is introduced by m 2 , p or D. This yields terms of the type 



d 4 xd A p 



tr 



-75 5£JV> If Nm $ Nm tyNp 1 ^Np 1 y/Np 1 + ■ 



+ o.t. 



(3.1) 



Here we have set Tj) R = <j) and p L =<j) + l^and ^represents the off diagonal (charged) part 
of p L as a matrix in ud space. We have kept terms with four W's and one derivative. Also 



we have introduced the quantity N = {p 2 



m 



2\-l 



. The dots in eq. (13. ip refer to further 



terms of the same type, while "o.t." refers to other terms which cannot have a contribution 
to the pattern 5Z(W + W-) 2 D. 

Next, we expand the ud labels using eq. (I2.2p for m and and eq. (I2.5P for 8Z. This 
produces 



6r 



d 4 xd 4 p 



tr 



-75 5ZN U i> \fr + CN d m d $ N d m d W~C~ l N u f> W + CN d j> VfTC^N^ + • 

(3.2) 



+o.t., 



where N U)d = (p 2 + m 2 u d ) 1 . 

At this point we can already factorize the trace between quantities which act only in 
generation space, namely, N u>d , m Ujd and C, and all the other quantities, which do not act 
on that space: 

. d^xd^p (\ 

oT 



( 9 _ ) . v -tr [N u CN 2 m 2 C- 1 N u CN d C- 1 N u }tr [75 5Z j> W + ? W~ f W + i> W~ i ] +• • • j +o.t . 

(3.3) 

It is a general rule that m u or m d can only appear raised to even powers and can be 
eliminated in favor of N u and iY^. 2 Eventually, all required momentum integrals and traces 



on 3 x 3 matrices in generation space can be cast in the form 

d 4 p 



16] 



a,b,c,d 



(27T) 



-{p 2 ) k ti[N a u CN b d C- 1 KCN d d C- 1 ] 



(3.4) 



where the exponents k, a, b, c, d are non negative integers. On the other hand, only the 
CP-odd contribution is of interest to us. This is the component antisymmetric under the 
exchange C — > C*, 



I a ,b,c,d — i I m I a ,b,c,d- 



(3.5) 



2 The label u or d does not change between two consecutive W's. This implies two consecutive Dl and so 
an even number of m's. That the presence of m u ,d can be obviated follows also from eq. (7.1) of 161 ] ■ 
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Due to cyclic and hermiticity properties of the trace and matrices involved, these integral 
satisfy 

Ia,b,c,d = ~Ic,b,a,d = ~K,d,c,b ■ (3-6) 

Such antisymmetry under exchange of the labels a and c, or b and d implies that many terms 
in eq. (13 .3p do not have a contribution to the CP-odd sector and this greatly alleviates the 
amount of subsequent computation. 

Performing an angular average over the momentum and taking the color and Dirac traces 
in eq. (13 .3p yields then, for CP-odd terms (N c = 3 is the number of colors) 

ST~ =N c j d 4 x (~/?, 1>2i2 5Z^W-d a W;W^Wt + •••)+ o.t. (3.7) 

It only remains to apply the derivative on all fields at its right. It can be checked that 
after this operation is carried out all terms so obtained cancel. So there is no term of the 
type (W + W-) 2 ZD in the CP-odd, P-odd sector of the Standard Model, in agreement with 



the alternative and more systematic calculation in 



16|. 



IV. DIMENSION EIGHT OPERATORS 



In this section we show that at eighth order in the derivative expansion there are non 
vanishing contributions in the CP-odd and P-odd sector of the Standard Model. Concretely 
we consider terms of the form (W + W~) 2 Z 3 D, with no gluons nor derivatives of the Higgs 
field, and apply the technique just described. In this case we select terms with n = 2 
in eq. (I2.10p and seek terms of the type 5Z(W + W~) 2 Z 2 D. The calculation is analogous 
to the one shown previously, except that now Z is not set to zero in p L , instead we use 
p L =(f) + tyf+ft. Then we keep terms with precisely four W's, two Z's and one derivative. 
After restoring the L, R labels and u, d labels, and carrying out the momentum integration, 
the trace in color, Dirac and generation space, and applying the derivative to the right, one 
obtains: 

6T- = N C J d 4 x(2lf X2A e^ Q p 5Z^Z x Z va W+WtW;W; + •••)+ o.t. (4.1) 

Here Z ua stands for the v derivative of Z a . We have eliminated the in favor of I\ hcd 

and have used identities involving 5^ u and e^ap to bring the expression to a canonical form. 
The expression contains 62 operators, each one weighted with various integrals I^cd- 

8 



It remains to find out the effective action from which the variation in eq. (14. ip derives. 
This serves also as a non trivial check of the computation. The method is just to propose 
all allowed independent terms of the form (W + W~) 2 Z 3 D with arbitrary coefficients, and 
take a first order variation with respect to Z to fix those coefficients. The Minkowski space 
result (see [if]] for further details in the conventions) is 



•l 
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;i2/ 1 - m 2 )z M zlw+w+w-w^ 



+ (4h + 23I 2 )Z^+W+WpW- + (6/i - 23I 2 )Z lx ZlW+W^W- 2 
+ (32/! + Ai 2 )Z^Z x Z a W^W^WpW; + (16/a - 38i 2 )Z fl Z x Z tT W+W+ W^W~ 
+ (16/! + 22I 2 )Z lx Z x Z (T W+W+W-W-p + (-20/ a - 15I 2 ) Z 2 x Z a W+W+W~W- p 
+10I 1 Z lx Z x Z ua W+ 2 WpW x - - 20I 2 Z fl Z x Z ua W+W+WpW- + c.c] (4.2) 

We have defined 1% = Ifi 2i — If 1 42 and I 2 = If 22S — I 21 3 2 , and "c.c" refers to complex 
conjugate; Z M is real, (W^)* = and I\ >2 are imaginary. The (underivated) Higgs field 
has been restored using that it scales as the mass dimension of l\ 2 . Also the derivative 



includes the field when it acts on the Ws [16j. Numerically, 3 

h,2= ( \\ 2 K 2 \ : «i = 0.226, k 2 = 0.456. (4.3) 

We can see that the values of these coefficients are considerably larger than simple esti- 
mates based on the Jarlskog determinant divided by the appropriate power of v . At sixth 
order the enhancement is driven by the small mass of the light quarks and so this can be 
considered as a kind of chiral enhancement. The possibility of such an effect was first pointed 



out in [ijj and confirmed in [15| 

The momentum integrals I abed are completely explicit but rather complicated homoge- 
neous functions of the fermion masses fl6j |. At sixth order these integrals are not continuous 
at rh u , m d , m s = 0, yet one can take the limit m n , rh d — > and subsequently m s — > and 



16]. 



this approximation gives a value fairly close to the exact one [16| . To discuss the situation at 
eighth order we will consider the simpler case of m , m t — > oo which is a quite good approx- 
imation for K\ 2 . At eighth order the momentum integrals are more ultraviolet convergent 
and also more infrared divergent than at sith order. Specifically, I 2 diverges as l/m 2 when 
fn u = frio- = 0, with k 2 = 1/2. The other integral, Ii, is more infrared divergent: in the same 



rn u ,md,m c ,m s ,mt,mb denote the quark masses. 
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limit K\ depends on the ratio fh u /fhd, varying continuously between —1/6 for fhd <C fh u to 
3/2 for fh u <C fhd . In eq. (14. 3[) we have used m„ = 2.55 MeV and fhd = 5.04 MeV. 

We have also considered CP- violating terms containing only W's and gluons. Such terms 
appear for the first time at eighth order since (at one loop) at least four W's are needed 
to violate CP and two G^ v are required to make a color singlet. In the P-odd sector one 
can write three independent operators, however we find that they have zero coupling in the 
Standard Model. On the other hand, in the P-even (C- violating) sector, there are also three 
operators of which one has zero coupling while the other two terms result in the following 
effective Lagrangian (in Minkowski space) 4 

£ eff (*) = -tph,2,2 (W^ 2 W-W~G; a Gl a - c.c.) . (4.4) 

Numerically, l\xii = iJ {^ir) 2 m 2 m 2 cl with n 3 = 3.76. This coefficient diverges loga- 
rithmically as fh u , fhd — > 0. Note that, at the order considered, the dimension four gluon 
condensate 5 does not induce a CP-violating interaction between the four W's. Such term 
vanishes identically, as it should be, since no CP-odd term can be written using just W's 
without derivatives or other fields. 

V. CONCLUSIONS 

We have shown that, to one-loop and at zero temperature, the leading P-violating CP- 
odd operators in the effective action of the Standard Model are of dimension eight. We 
have computed explictly the couplings for the operators of the form Z 3 (W + W~) 2 plus one 
covariant derivative, eq. (14. 2p . These operators come from Feynman graphs with one quark- 
loop, four W legs, three Z legs and possibly one B ~ Z + 7 leg. In principle, dimension 
six operators could develop beyond one-loop or at finite temperature due to the breaking of 
Lorentz invariance. Purely gluonic operators of dimension eight have also been computed, 
eq. (I4.4p . and they are C-odd and P-even. Remarkably the corresponding coupling constants 
we find are not vanishingly small, rather they have a natural scale related to intermediate 
mass quarks times the Jarlskog invariant. All formulas derived for quarks extend directly 
to massive Dirac leptons. This implies that even if the neutrino masses are small their 

4 The gluon field strength tensor has been normalized according to [D^D,,] = £(A a /2)G" . 

5 In the presence of the gluon condensate G^G^n = ^(5 Ma £ K/ g — ^/9^c)(( ( ^Acr) 2 ) + fluctuations. 



10 



contribution to the CP-violating couplings needs not be small, due to infrared sensitivity in 
the momentum integrals on which the couplings depend. As a consequence, such couplings 
will be strongly dependent on the mass ratios between neutrinos of the different generations. 
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